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MEMORANDUM 5-5-59W 

EFFECT OF CONTRACTION ON TURBULENCE AND TEMPERATURE 
FLUCTUATIONS GENERATED BY A WARM GRID 
By Robert R. Mills, Jr., and Stanley Corrsin 

SUMMARY 


Hot-wire anemometer measurements were made of several statistical 
properties of approximately homogeneous and isotropic fields of turbu- 
lence and temperature fluctuations generated by a warm grid in a uniform 
airstream sent through a k-to-1 contraction. These measurements were 
made both in the contraction and in the axi symmetric domain farther 
downstream. 

In addition to confirming the well-known turbulence anisotropy 
induced by strain, the data show effects on the skewnesses of both lon- 
gitudinal velocity fluctuation (which has zero skewness in isotropic 
turbulence) and its derivative. 

The. concomitant anisotropy in the temperature field accelerates the 
decay of temperature fluctuations . 


INTRODUCTION 


Significant measurements of the effects of gross strain upon homo- 
geneous turbulence have been made by Townsend (ref. 1) and Uberoi 
(refs. 2 and 3)- Townsend ran grid-produced turbulence through a nearly 
constant-area distortion in which the channel section went from a L-to-1 
rectangle to a 1-to-k rectangle. The mean velocity was nearly constant. 
Uberoi ran grid-produced turbulence through a square channel contraction, 
which was essentially an axially symmetric distortion. He tested with 
contraction ratios of b to 1, 9 to 1, and l6 to 1. 

They both made measurements of the three components of turbulent 
kinetic energy, the mean-square derivatives (giving the values of dissi- 
pation scales), as well as of the pertinent turbulent energy spectra 
before and after distortion. 
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Their principal results may be qualitatively summarized as follows: 

(a) Distortion of roughly isotropic turbulence increases the veloc- 
ities associated with the vorticity components which amplify by vortex 
line stretching. This phenomenon was predicted and estimated by Prandtl 
(ref. b ) . 

(b) The gross strain destroys isotropy not only of the energy -bearing 
eddies but also of the turbulent fine structure. There was departure from 
"local isotropy" in the distorting part of the channels. 

(c) The longitudinal and lateral spectra are deformed somewhat as 
predicted by the calculations of Ribner and Tucker (ref. 5) and Batchelor 
and Proudman ( ref . 6 ) . 

(d) Neither arrangement corresponded to a "rapid" distortion. There 
was appreciable transfer of energy among turbulence components during the 
straining process. 

(e) In a straight channel following tha distorting section, the 
return to local isotropy is fairly rapid, bat the return to general 
isotropy is relatively slow. 

Direct theoretical work on this problem began with Prandtl who used 
models consisting of steady, invi sc id -stream tubes and vortex tubes 
(ref. 4). This work was followed by the ce Llular vortex analysis of 
Taylor (ref. 7)> which is the basis for the more recent work. Ribner 
and Tucker (ref. 5) and Batchelor and Proudman (ref. 6) independently 
applied Taylor's model more realistically by simple Fourier synthesis. 

This most recent theory assumes (l) instantaneous distortion, so there 
is no decay or intercomponent transfer during distortion, and (2) homo- 
geneous distortion. Ribner and Tucker appended a semiempirical correc- 
tion to include decay when the distortion i; not instantaneous. Using 
the same model, Acharya recently calculated some details of the case in 
which the unstrained field is axisymmetric 'ref. 8). 

The present investigation was undertaken to confirm and extend the 
Townsend -Uberoi findings on turbulence dynamics. In addition, a roughly 
isotropic (scalar) field of small temperature fluctuations was added to 
the turbulence so that the effect of gross distortion upon homogeneous 
turbulent mixing could be studied. Explicit predictions on the relative 
structure of scalar and velocity fields have been given chiefly in the 
isotropic case (refs. 9 to 11), although some exploratory analysis has 
also been directed at the axisymmetric. case (refs. 12 and 13). Experi- 
mental study of this homogeneous mixing problem has been confined to the 
nearly isotropic fields behind a warm grid 'refs, ll and 15) . 
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It has long been known that in the absence of shear stresses aniso- 
tropic turbulence will tend to become isotropic at sufficiently large 
Reynolds numbers. For example, the turbulence generated by a regular 
grid in a uniform stream reaches a fairly homogeneous but somewhat axi- 
symmetric state by perhaps 30 mesh lengths downstream (refs. 16 and 17). 
This tends to become isotropic farther downstream. Townsend and Uberoi 
found this tendency in their postdistortion fields. Remarkably, the 
theoretical studies of homogeneous axisymmetric turbulence (refs. 18 to 
20 and 12) have still produced no prediction of the rate at which this 
occurs. Therefore, it seemed desirable to make some detailed measure- 
ments at the downstream end of the channel for comparison with the 
immediate postcontraction turbulence. 

This work was carried out at the Johns Hopkins University under the 
sponsorship and with the financial support of the National Advisory 
Committee for Aeronautics. The authors would like to thank Dr. L. S. G. 
Kovasznay, Dr. R. Betchov, Dr. F. H. Clauser, Mr. L. T. Miller, and 
Dr . D. W. Dunn of the Aeronautics Department of the Johns’ Hopkins 
University for their helpful interest and Mr. S. Bhaduri for his care- 
ful preparation of the figures. 


SYMBOLS 


a 

c 


E 


ei,e 2 


hot-wire sensitivity 


mean velocity ratio (hence 
(when used as a variable,, 


streamvise distortion) 
c(x) = 


y 



one -dimensional energy spectrum (subscript indicates field 
variable) 


voltage fluctuations across two hot wires 


f,g principal two-point velocity correlations in isotropic 

turbulence 


K 

k 

L 


two-point velocity triple -correlation coefficient of u, 

u 1 2u 2/ /ui 2 u 2 ' 

wave number in x-direction 

integral scale of random field (subscript indicates field 


variable) , 


0 


(Correlation function)dr 



k 


M 

m 


n 


grid mesh size 

two -point temperature double -correlation coefficient, 
0-^02 ? 

two-point temperature triple -correlation coefficient, 

9 l 20 2/ 9 l 20 2 , 


m,n indices 

p static pressure 

q two-point velocity double-correlation-coefficient functions 

q resultant fluctuation velocity, v/u^ + v^ + w 2 

r scalar distance between two observation points 


SjS^Sq skewness factors for 


du 

bx.' 


u, and 0, respectively 


t time 

U,V,W mean velocity components along x, y, and z axes 

u, v, w velocity fluctuation components along x, y, and z axes, 

respectively 


x,y,z Cartesian coordinates 

x • > x • i X t, Cartesian coordinates, tensor notation 

a angle between tunnel axis and lire joining deflection points 

^ in section downstream of contrection 


p angle measured around tunnel axis 

y thermal diffusivity 

q coordinate difference in y-direcl ion 

© mean temperature rise across grid 

0 temperature fluctuation 

\ dissipative scale of random field (subscript indicates 

field variable) 
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v kinematic viscosity 

I coordinate difference in x-direction 


p 

fluid density 


n 

average value 


( r 

root -mean-square value y ^ 


Subscript 

s : 


ax 

axisymmetric 



refer to f and g correlations 

9 respectively 

j 

tensor notation referring to x-, 
respectively 

y-, and z-directions, 

iso 

value for isotropic (undistorted) 

field at equal time 

l ,m,n 

indices 


0 

reference value 



u,v,6 refer to u, v, and 6, respectively 
1,2 upstream and downstream of contraction, respectively 

•I parallel to tunnel axis 

■1 perpendicular to tunnel axis 

EXPERIMENTAL ARRANGEMENT 

The basic wind tunnel was the open-return unit sketched in refer- 
ence 9- The only addition was a secondary contraction with 4-to-l veloc- 
ity ratio starting 19^ inches downstream of the grid (fig- 1). After the 

8 

contraction was an almost uniform channel, slightly divergent to maintain 
constant static pressure (hence constant velocity) in spite of boundary- 
layer growth. This channel permitted some study of the history of the 
axisymmetric turbulence emerging from the contraction. 
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The measured mean-velocity distribution of the tunnel is shown in 
figure 2. The free-stream velocity fluctuation levels upstream of the 

grid were — « 0.08 percent and i- = ^ 0.1 percent. An apprecia- 

u u u 

ble part of this was pulsation resulting from the unbalance inevitable 
in the low-quality centrifugal fan used to move the air. The relative 
excess energy in lateral fluctuations is, of course, characteristic of 
wind tunnels with contractions and is consistent with the behavior of 
the contracted grid turbulence which is the primary subject of this 
study . 

The biplane, square-mesh, fluctuation-generating grid was of 
l/4-inch -diameter, round metal rods (Calrod heating units) set 1 inch 
on centers. It was heated with 60-cycle, 220-volt, 5-phase alternating 
current, and had sufficient thermal inertia that no 120-cycle periodic 
temperature fluctuations could be detected downstream. 

A principal difficulty to be overcome was a "built-in" vertical 
mean-temperature gradient characteristic of the wind-tunnel-plus-room 
system. Individual current controls on the horizontal rods permitted 
reducing this to a difference of 0.2° C over the height of the tunnel 
working section. Figure 8 of reference 15 shows this, as well as a 
typical horizontal distribution of mean temperature, which was uniform 
to within ±0.05° C. The mean-temperature rise through the grid 0 was 
5° C at a wind speed of 14 feet per second. 

Both velocity and temperature fluctuations were reasonably homoge- 
neous over the tunnel cross section at some distance behind the grid. 
Topical traverses of root -mean -square temperature 9' and velocity u' 
fluctuations before and after contraction are given in figures 5( a ) and 
5(b). The two sets of points in figure 5(t) were run at different times 

The negligibility of the effect of air-density fluctuations upon 
the turbulence dynamics was confirmed by tie equality of turbulence 
levels and velocity correlation functions with and without heating. 


MEASURING EQUXPMINT 


Velocity and temperature fluctuations were both measured with a 
hot-wire anemometer (ref. 21). The basic equipment is described in 
reference 22. The auxiliary circuitry for squaring, cubing, etc. is 
described in reference 15* 

Because of the extremely small temperature fluctuations to be 
detected, unusually high-resistance hot-wires (between 75 and- 150 ohms) 
were used. These were of 0.00005-inch-diameter platinum (Wollaston) 
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about 0.1 inch long. No length corrections have been applied to the 
data, but the high wire resistance necessitated correction for current 
fluctuations in this nominally constant current equipment. 


EXPERIMENTAL PROCEDURES 


, .e experimental procedures used here were largely the same as those 
escribed in reference 15 . The only additional variables were the lat- 
eral velocity fluctuations v and w. These were detected in standard 
as ion by the difference voltage of a (directionally sensitive) x-meter. 

Power spectra of u, v, and 9 were determined with the aid of a 
Hewlett-Packard 300-A wave analyzer. All double correlations as well 
as the isotropic precontraction triple correlations were obtained as 
in reference 15, but the postcontraction triple correlations of u 

required a bit more elaborate procedure in regions where v? 0. The 
extra effect can be schematically shown by assuming linearized and iden- 
tical wire response, compensated for thermal lag, so that the two output 
voltage fluctuations due to velocity are 


e l = -au(x,y,z,t) 
e 2 = -au(x+| ,y+t], z, t ) 


Then 



In true isotropic or axi symmetric turbulence , u^ must be zero^ so 
is also zero within the approximation of linear wire response. 

In practice, of course, equation (2) was not used directly but was 
corrected for (a) unequal sensitivities of the two wires and (b) non- 
linearity of wire response, as outlined in the appendix of reference 15. 

An attempt was made to determine the triple correlations 

0 2 (x,y,z,t)0(x+^,y+r 1 ,z,t) and u(x,y,z,t)0(x,y,z,t)0(x+|,y+ri,z,t), but 
the combination of small thermal signals and complex triple-correlation 
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computation procedure yielded excessive scatter in the data, especially 
for the latter quantity. 


EXPERIMENTAL RESULTS 


Figure 4 demonstrates that the velocity and temperature fields are 
both axisymmetric after the contraction. This shows double-correlation 
values measured with two wires separated by a fixed radius in the same 
yz plane (a-2 = 90°). One wire was fixed on the tunnel axis, and the 
other was rotated to a series of positions on a circle about the axis. 
The value of r is the radius of this circle. 

Figure 5 ( a ) gives local turbulence levels along the tunnel axis. 

It was found that w 1 ~ v ! , so only v f data were taken extensively. 

A more dramatic indication of the amplification of lateral fluctuations 
is given by numbers proportional to the root -mean- square fluctuations 
themselves, as in figure 5("b) • 

For comparing the histories of strained and unstrained turbulence, 
time rather than distance is obviously the appropriate reference. The 
conversion curve from x/M to a dimensionless time tU]_/M is given in 

figure 6. Figure 7(a) shows the first such comparison, the square root 

of the total turbulent energy (^q 1 = \/u^ + v^ + w^ ). The unstrained 
results here, as in several other figures, are taken from reference 15- 
These results are sometimes called isotropic. The data on root-mean- 
square temperature fluctuation (fig. 7(b)) are contrasted with the 
unstrained case in figure 7(c). Evidently the distortion hastens 
decay. 

A class of important statistical poin: functions of turbulence is 
the ’’skewness factor” of any variable. This is a dimensionless measure 
of the third moment of the probability density and is central in that it 
may be a measure of nonlinear effects. In isotropic turbulence, the 
velocity fluctuations themselves must have zero skewness because of the 
required invariance to coordinate reflection; for example, 



This condition is well met by grid-generated turbulence. However, isot- 
ropy does not require that velocity-component derivatives in the direc- 
tion of the components du/dx, 8v/8y, and bw/bz have zero skewness. 
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One of the quantitative challenges to any isotropic theory is the pre- 
diction of the generally measured value of about -0.4 for qua n tities 
such as 


S = 



8(a) and 8(b) show the considerable effect of the contrac- 
tion on these two skewness factors. Figure 8(c) summarizes the data on 
a time scale. 

The small discrepancy between S in the unstrained field just 
ahead of contraction and S at the same station without later contrac- 
tion (ref. 15 ) is probably a chance inconsistency between sets of data 
taken some weeks apart and at slightly different locations in the cross 
section of the duct. 

The return to isotropic values by both Su and S is evidently 
more rapid than the return of v'/u 

A condition of isotropy in a scalar field like 0 does not require 
that 0 have statistical symmetry. Nevertheless, it was found in the 
unstrained case (ref. 15 ) that Sq approaches zero rapidly. Figures 8(d) 
and 8(e) show that although the strained 0 approaches zero skewness 
later in space, it actually goes down faster in time. 

The temperature-gradient fluctuation, being a vector, must be sta- 
tistically symmetrical in an isotropic field, and an attempt was made to 
measure the skewness of de/dx along the axis of the duct. Unfortu- 
nately, differentiation with respect to time ( — — ** decreased the 

\dx u btj 

signal-to-noise ratio too much for reliable measurement. 

The effect of contraction on the energy spectra of longitudinal and 
lateral velocity fluctuations is at least qualitatively detectable in 
figures 9(a) and 9(^0 • The abscissae of the postcontraction curves are 

U 2 

— k instead of k so that the sinple geometrical distortion effect is 
Ui 

removed. Relatively, then, it is clear that the strain has favored the 
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small structure in u and the large structure in v. This is consistent 
with previous work, both experimental and theoretical. 

The temperature spectra, similarly compared (fig. 9(c)), show less 
difference. This means that the contraction has acted pr im arily as a 
geometrical distortion of the scalar field. 

Double- and triple-correlation functions downstream of the contrac- 
tion at x/M = 4-5 and 73 (fig- 10) provide a fairly detailed picture of 
the nature of the distorted fields and their subsequent development. The 
precontraction correlations at x/M = 17 are given in figures 17(b), 
19(h), 2l(a), and 22(a) of reference 15- 

The double correlations of u and 0 are each given for three dif- 
ferent values of oq, the angle between the tunnel axis and the line 
joining the detection points. The v-correlations were measured at only 
two angles. Comparison of figures 10(a) and 10(b) showfe a decided trend 
toward isotropy in the straight duct following the contraction. 

Typical two -point velocity triple-correlation functions in the post- 
contraction turbulence (figs. 10(g) and 10(h)) also show a trend toward 
isotropy. Inadequate signal-to-noise ratio Led to excessive scatter in 
the attempts to measure the temperature triple correlations. Still, the 
general character for three different directions after distortion at 
x/M = 45 can be gleaned from figure 10(i). The corresponding precon- 
traction data (nearly isotropic, at x/M = 17) are in figure 25 of 
reference 15- 

"integral scales" are defined from velocity double-correlation 
functions. At a given (nominal) value of x: 


III 

f Q||( r ^?) dr 

111 

l 

hi" 

O 00 

/ Q±( r ^ a : > ) dr 

J 0 


(3) 


In figure 11(a) these are 

scales. 


compared with the isotropic integral 


Lf = 



f(r) dr 


(*) 


L, 


0 


g(r) dr 
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where, in isotropic turbulence, 

Q||(r,0°) = Q ± (r,90°) = f(r) 

( 5 ) 

Qj_(r,0°) = Q||(r,90°) = g(r) 

J 

Therefore, in figure 11(a), L,|(5°) and L x (90°) have departed from 

the L f curve. There is not enough information in the figure to tell 

even roughly when these integral scales will reach their isotropic rela- 
tive values. 


The same is true of the thermal integral scales (fig. 11 (b)). In 
this case, however, the rough spectral similarity maintained through the 
contraction (fig. 9(c)) suggests that one look for a corresponding pro- 
portional distortion of the values of L e : at ijjk = 27, ^ - 


and 


L 0 (9O°) 


L0 


iso 


7 iso 


0.6. For simple instantaneous affine distortion these 


values would be b and 0.5, respectively. 


Some of the dissipative scales have been obtained graphically from 
double -correlation functions. They are the abscissa intercepts of the 
vertex -osculating parabolas : 


A|, 2 (a2) 


= 1 
~ ~ 2 


dr' 


2 Q l|( r ' a 2) 


r=0 


\L 2 (a£) 


s _ 1 
2 


is ^( r '°2> 


r=0 


2 

Aq (op) 


dr 


2 m ( r > a 2) 


J r=0 


(6) 


Hie values of A ( | and A, are compared with the undistorted dis- 
sipative scales in figure 12 (a) . These isotropic values are taken from 
reference 15: 
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V 


1_ = l(d‘ f 

2 


2 W 2 /r=0 


17 = 5 WU, 


r=0 


(7) 


Here A||(5°) and. A _^( 90° ) are departures from Af ; Aj|(90°) and Aj_(5°) 
are departures from A„- As with the integral scales, not all of these 

O 

quantities appear to return rapidly to their isotropic relative values. 
However, the more extensive data on mean-square velocity derivatives 
(also related to the values of A) taken by Townsend and Uberoi do show 
a strong return to isotropy. The evidence In figure 12(a) must be 
regarded as simply inadequate to permit any independent conclusion. 


tUi / , 

The values of 'Kq at -jjp = 27 (fig. 12(b)) are moderately con- 
istent with interpretation in terms of simple geometrical distortion 

M5°) 


through the contraction: 
tU]_ 


^®iso 

tu 


k.l 


Aq(90°) 

and — 1 « 0 . 7 . 


*9 


In the 


iso 


interval from _ , 

M M 

toward isotropy in the Aq values. 


UU 1 1 

= 27 to = there is no appreciable return 


Figure 13 is a rough estimate, from the data, of the relative decay 
rates of temperature and velocity fluctuations. In the unstrained case 
(dashed line, from ref. 15 ) this quantity appears to approach an asymptote 
of 0.6. 


TURBULENCE LEVELS 


The behavior of v’/u* in a contraction is qualitatively given by 
Prandtl's theory of differential vortex stretching (ref. 4). For a quan- 
titative theoretical prediction one can turn to the Ribner-Tucker analy- 
sis (ref. 3). Values obtained by this analysis, including a semiempirical 
estimate of energy decay from isotropic data, compare with those in fig- 
ure 3(b) as follows: 
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Ratio 

Ribner-Tucker 

method 

Measured 

v/v 

0.4 

0.7 

v 2 ' /vx ' 

1-5 

1.4 

( v 2 ' / u 2 ' ) 
(V/V) 

3-8 

2.0 


Only the v' prediction appears to be satisfactory, but even this 
conclusion must be tempered by the likelihood that there is appreciable 
transfer of energy from v' and w' to u' during the straining proc- 
ess. Such an effect could account for the particularly low u 2 ' esti- 
mated by the Ribner-Tucker method. On the basis of total turbulent 
energy, the agreement is good: 


Ratio 

Ribner-Tucker 

method 

Measured 

<12 7*1 ' 

1.2 

1.2 


It was concluded that the gross strain rate inflicted on the turbu- 
lence field was far from rapid, a conclusion also reached by Townsend 
and by Uberoi for their experiments. Both the energy decay and the inter- 
component transfer were appreciable during the 0.033-second-contraction 
transit time. A characteristic decay time is >\g/u', which was about 

0.025 second just entering the contraction. For a rapid strain, the 
transit time would have to be much less than this . 

There is still insufficient axisymmetric turbulence theory to define 
a time characterizing its rate of approach to isotropy, for example, an 
intercomponent -energy transfer time. If, however, this is an inertial 
effect (acting through the pressure-velocity correlations), L g /u' might 
be a plausible choice. This was 0.04-4 second entering the contraction. 

The foregoing paragraphs raise the question of whether it is feasi- 
ble to produce rapid contraction of turbulence in the laboratory. For 
example, the contraction nozzle used in the present tests probably could 
not be shortened by more than a factor of 2 without inducing boundary- 
layer separation immediately upstream or downstream. But even halving 
the transit time would not make it negligible compared with \?/u' or 
Lg/u 1 . 
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Foregoing the assumption of instantaneous distortion, the down- 
stream development of the turbulence can be displayed by means of 
Reynolds' turbulent -energy equations. 

For the three component energies, 



A first approximation, suitable for the present case, follows from 
the assumption that turbulent mean values are nearly constant in the 

|du 2 l 


yz plane, for example, that 


b\J 

« 

dU 

dy 


dx 


du^ 

dy 


« 


dx 


This also involves 


, although there is no such inequality for V and W. Then 


1 tt du 2 _ ,,2 dU 1 i — o- 

5 U 5T'- U ^-55T'p u S +vu ^ u 


( 12 ) 


t“tr 


-V 2 ^ - + vv^v 

by 2 P c)y 


(13) 




-W 


2 M - 1 - lw^ + wV 2 * 

£z 2 dx P + 


(14) 
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In "this approximation , it is also assumed that uv, vw, and wu are 
negligible . Adding to get the total turbulent kinetic energy and usi ng 

axial symmetry (e.g., v 2 = w 2 ) as well as mean -flow continuity. 



(15) 


A still rougher approximation follows from restricting to very small 
turbulence levels, so that one can expect inequalities like 


du^ 

« u I du2 

dx 

j dx 


(16) 


If p' is of the order pu 2 , the static -pressure term can also be 
neglected : 


1 u » _ u 2 &2 + vu ^ u 

2 Sx c* 


(17) 


Since v 2 = w 2 , add equations (lj) and (l4) and eliminate V and W in 
terms of U : 


S “ P - £ ^ 

^ OX 2 dx 


(18) 


Comparison of the first terms on the right-hand sides of equations (17) 
and (18) shows clearly that values of > 0 tend to decrease u 2 and 
to increase v 2 and w 2 . 

The q2 equation simplifies to 



(19) 


Equations (17) and (18) cannot be used further without an estimate 
of the viscous terms. First, it is convenient to transform 


uV^u = i V^u^ - 
2 


chi chi 


( 20 ) 
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with corresponding forms for v, w, and u j . Even in turbulent shear 

flows the first of these terms is usually negligible, and the second is 
a good approximation to the dissipation. Then equations ( 17 ) and (l 8 ) 
are 


1 T T Ou^ 

2 U cbT 


_ u 2 Ou v du Ou 

u v ^r^r 


( 21 ) 


1 u M « i v2 ^ - v 8v_ 0v_ 


Ox 


Ox 


C»Xj dxj 


(22) 


An estimate, very rough indeed, could now be made by assuming simple 
geometrical distortion to relate these viscous terms to those in unstrained 
turbulence at equal times. It is already well established, however, that 
the behavior of the velocity field under strain involves dynamics as well 
as kinematics. It seems more reasonable, therefore, to save the technique 
outlined for a more appropriate problem, tiat of the strained temperature 
field. 


TEMPERATURE FLUCTUATIONS 


Without mean-temperature gradients, fluid strain has no effect on 
the probability density of the 9 -field. Hence it has no instantaneous 

effect on moments such as 9 ^. The strain obviously changes the 

temperature-gradient field, however, and hence the rate of decay of 9 ^. 

In fact, it is just this kind of effect which enables turbulent 
motion (through its random strain field) to promote rapid molecular- 
scale mixing. 

The crudest estimate of mean-strain-aocelerated 9^ decay (fig. 7(c)) 
can be made by assuming statistically homogeneous and uniform mean strain. 
Then the decay equation is (ref. 10) 


d9^ _ 89 d(i 

dt Oxi 


( 23 ) 


Suppose the field is instantaneously strained by a factor c [= — j in 

V Ul / 

1 / Vc in the y- and z-directions . With an isotropic 


the x-direction and 
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field before distortion, 



postdistortion derivatives are 



and the immediate 



(24) 


The decay rate immediately after distortion is then 


d9^ 

dt 



(25) 


Assuming the unstrained history 0 ± 2 (-t) is known, the relative decay 
rate is initially 


d9 2 

de. 2 

1 



(26) 


Rather than attempt to apply this instantaneous -distortion analysis to 
the experiments, it seems worthwhile to consider a more general approach. 

A more reasonable estimate requires (a) using the general 02 equa- 
tion without assuming homogeneity (ref. 23), and (b) including c = c(x) 
where c(x) = U(x)/Up. With constant mean temperature. 


5 Ui if - - § + ^ 


With x the direction of largest gradients, 

I " 57 “ - 1 ^ 


2 8x 


(27) 


( 28 ) 
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and for very small turbulence level, the triple correlation term is 
neglected. Making use of a transformation like equation (20) and 
neglecting the second derivative term, 

U « _2y £i_ <*L_ (29) 

Ox dxi 8xi 


Finally, ignoring the differences between the diffusion-induced growth 
rate of in strained and unstrained flows. 



where c(x) = u (2il and Aq . is the dissipative scale for the unstrained 

Ul 1 

case at the same time . 

With time t as the independent variable, 553 ■— and equa- 

tion (29) becomes 


d£ „ _h 7 JL(. i- + 2c 
dt Tvq^Xc 2 


With given values of c(t) and (t) this could be numerically inte- 

grated, but the latter can be eliminated by seeking only the ratio of 0^ 
in the strained and isotropic cases. 

In terms of time, equation (29) can be written for the isotropic 

case as 


-127 




dGi 

dt 


( 32 ) 
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From equations (31) and. (32), 



Numerical integration of equation (34) from x/M = 15 (where it 
starts to deviate from the unstrained curve) to x/M = 35 gives as theo- 
retical estimate for the difference 


0i' - 0' 


9 ' 


0.15 


But the experimental value, from figure 8(e), is 


9 


i 



0.35 


which is more than twice as large. The weakest link in the analysis is 

probably equation (30) . Inclusion of the effect of increased q2 (which 
reduces dissipative scales by increasing turbulence Reynolds number) only 
increases the theoretical estimate to 0. 19- 


SKEWNESS AND INHOMOGENEITY 


For turbulence dynamics, the principal results are the effects of 


contraction on the skewnesses of u and 


8u 

5F 


i • The generation 


U 


of S^ during the straining process, when its value was negligible 

upstream, is evidence that the strain is appreciably inhomogeneous. The 
statistical symmetry of a homogeneous strain would preclude the produc- 
tion of u3 unless it were present at the start. 
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A good theoretical estimate of u3(x) is made difficult to obtain 
by the lack of knowledge of pressure-field statistics. The equation for 

u3 is derived by multiplying the u-equation by u 2 and then averaging: 


lu^ + iv^+iw|^ + u5|J + u 2 v ^ + u 2 w |W + i^ + I c^3v + 

3 dx 3 oy 3 oz dx qy dz 3 ox 5 dy 


1 du3w U P du 2 U P duv U P duw _ - i — + V35BS (55) 

3 dz dx dy dz P Ox 


For a first approximation, it may be reasonable to assume 


j 

j 


l-l 


v| = 

1 w 

« 

M 


du3 


du3 

« 

di.3 

dy 


dz 


| by 

U 2 V 

= 

U^w j 

« 

3 


uv = uw 


« u 2 


(36) 


Then equation (35) reduces to 


iu^~-u3^-i^: + i - I u 2 ^ + vu 2 ^u (37) 

3dx dx 3 dx 2 <^x pdx 


Restriction to a low turbulence level does not permit neglecting the u^ 

and (u 2 ) terms relative to the u3 U teims because u3 is initially 

zero while u 14 ' and u 2 are not. On the other hand, u 2 ^ is zero in 

dx 

isotropic turbulence. If it is assumed that p 1 = O^pu 2 ^, restriction 
to a small turbulence level permits neglecting the pressure -velocity 
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term. If, furthermore, 7 + is estimated roughly in terms of u 2 by 
taking the Gaussian equality u 4 13 3 (u^) , equation (36) becomes 


in du3 ^ 3 l (u^ ) o, o 

3 u sr u 5 - 2 "V- + ^ 


(38) 


Some qualitative confirmation of this form can be obtained by using 
it to predict the sign of u3 in the early stages of the straining proc- 
ess. If u3 is identically zero at the origin of x and the effect of 
viscous force is ignored, equation (38) predicts u 3 > 0 for small values 

of x. In the present experiments it turned out that u 3 > 0 for all 
values of x. 

Doubtless some knowledge of the viscous term in equation (38), 


u^vu = — — hZ _ 8u 8u^ 


3 dx-^dxi 


dxj 8xj 


8u 8u^ 
dXj 8xj 


(39) 


would be required for an estimate of u 3 ( x ). On the other hand, there 

u3 


(u')- 


V 


may be less effect of viscosity on the dimensionless ratio 
Furthermore, it is possible that u 2 in equation (36) may still be 

OX 

negligible for low enough turbulence levels: with distortion it departs 

from zero, but so does the term u3 ^U. 

8x 


By chain differentiation. 


8Su _ 8 j 

u 3 

1 du 3 

3 $u du^ 

8x 8x 

/— \ 3 /2 

/— \ 3 / 2 8x 

^ ~ dx 


(u 2 ) 

M 



( 40 ) 
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Substituting from equations (17) and (38), less viscous terms, this 
becomes 


« 6 iC ^ (4i) 

dx ^ dx 


Furthermore, u' can be expressed in terms of U by integrating the 
inviscid form of equation (17) 


U '( X ) « Jit-up (42) 

U(x) 

so that 

65. u , 1 4- a ( 45 ) 

dx x U' “ 

Therefore 

s u ( x )» 5 VLii\ (44) 

Ui \ u (x)/ 

Entering the contraction, — 25 0.025- At the end, U/U^ = 4, so equa- 

Ul 

tion (43) predicts that S u * 0.07 at the end; the measured value is 
about 0.075* 

It is clear that this good agreement must be in part fortuitous. 

For example, in the straight section following the contraction equa- 
tion (45) predicts a constant value of £> a 55:5 0.07- Yet, in fact, 

returns rather rapidly to zero. It has long been known that the inter- 
component energy transfer (accompanying a trend toward isotropy in homo- 
geneous turbulence) depends upon pressure -velocity correlations which 
were neglected in this analysis. 
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CONCLUDING REMARKS 


Measurements of the turbulence levels of grid-produced turbulence 
passing through a 4-to-l contraction and a downstream straight duct agree 
with the results of Uberoi. They confirm in general the conclusions of 
both Townsend and Uberoi that reasonably smooth channel distortions can- 
not be regarded as sudden in terms of ordinary grid turbulence structure. 

In addition, the relatively strong effects of this contraction on 
skewnesses of both the longitudinal velocity fluctuation and its deriva- 
tive indicate that the contraction cannot even be regarded as homogeneous. 

Establishment of this contraction as a rather slow and inhomogeneous 
distortion of this turbulence reduces the incentive for carrying out a 
detailed comparison of the correlation functions with predictions of a 
Ribner-Tucker-Batchelor-Proudman type of analysis. Nevertheless, the 
homogeneous, instantaneous approach accompanied by correction for dissi- 
pation does seem to give a satisfactory estimate for total turbulent- 
energy change . 

The effect of a contraction in increasing the root -mean- square veloc- 
ity fluctuation components v ! and w' relative to u' suggests that 
nearly isotropic turbulence behind a grid (where v 1 and w’ are less 
than u 1 ) may be obtained by a suitably small amount of contraction prior 
to the straight duct. (Uberoi reported in a private communication that 
he has tried this and found a strong tendency for the turbulence in the 
straight section to return to its axisymmetric state.) 

The accelerated temperature-fluctuation decay due to strain is qual- 
itatively understandable. The failure of the present (noninstantaneous) 
simple geometrical theory certainly indicates the need for a more detailed 
look at the problem. 


Mechanical Engineering Department, 

The Johns Hopkins University, 

Baltimore, Md., March 23, 1958* 
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Figure 2.- Measured mean-velocity distribution of wind tunnel 








(a) Before contraction; x/ M = I 7 . 0 . 


Figure 3*- Topical traverses of root -mean -square temperature 9 

velocity u r fluctuations . 











Figure 5*- Longitudinal and lateral turbulence decay. = 7^20. 
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Figure 7*- Comparison of histories of strained and unstrained turbulence. 







(b) Root-mean-square temperature fluctuation for strained turbulence 
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(c) Curve from figure 7(b) contrasted with curve for unstrained turbulence plotted against 

dimensionless time. 8 = 5° C • 




Figure 8.- Effect of contraction on skewness factors. Different points indicate data taken 

at different times. 
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(c) Skewness factors for du/dx and u plotted against time 
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Figure 8.- Concluded. 
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(a) Spectrum of u^. 

Figure 9 .- Effect of contraction on energy spectra of velocity 

U2 , 

fluctuations, c = = 4. 
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(c) Spec 


Figure 9 
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Concluded. 
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Figure 10.- Double- and triple-correlation functions downstream of contraction 




k6 



Figure 10 . - Continued . 
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igure 10.- Continued. 
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Figure 10.- Continued. 
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Figure 10.- Continued. 
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Figure 11.- Comparison of integral scales of random fields. Curves taken from reference 15* 











(b) Thermal dissipative scale 






